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Abstract. A fundamental result of free probability theory due to Voiculescu 
and subsequently refined by many authors states that conjugation by indepen- 
dent Haar-distributed random unitary matrices delivers asymptotic freeness. 
In this paper we exhibit many other systems of random unitary matrices that, 
when used for conjugation, lead to freeness. We do so by first proving a general 
result asserting "asymptotic liberation" under quite mild conditions, and then 
we explain how to specialize these general results in a striking way by exploit- 
ing Hadamard matrices. In particular, we recover and generalize results of the 
second-named author concerning the limiting distribution of singular values of 
a randomly chosen submatrix of a discrete Fourier transform matrix. 



1. Introduction 



Of the results of Voiculescu [21] [22] providing the foundations for free probability 
theory, arguably the simplest and most familiar is the following. Let A^-^^^ and B^-^^^ 
be deterministic N-hy-N hermitian matrices with singular values bounded indepen- 
dently of TV and having empirical distributions of eigenvalues tending in moments 
to hmits and fiB, respectively. Let C/^^^ be an N-hy-N Haar-distributed ran- 
dom unitary matrix. Then the empirical distribution of eigenvalues of the sum 
yl(^) _|- j^{N)jj(N)* ^gQ(jg moments to the free additive convolution Hms- 
The question addressed here, very roughly speaking, is this: how much less random 
can we make J/^^-' and still get free additive convolution in the limit? More generally, 
we ask: what sorts of random unitary matrices J/^^-' make and [/^^-'S'^-'C/''^'* 
asymptotically free? That the randomness could be reduced a long way was sug- 
gested to us by the main result of [Hj concerning the limiting distribution of singu- 
lar values of a randomly chosen submatrix of a discrete Fourier transform matrix. 
When the latter result was discovered it seemed related to free probability theory 
only by coincidence; here (see Corollary 13.61 below) we show that the result can 
be well-understood as a special case of a much more general phenomenon in free 
probability. 

In a bit more detail, our explanation of the results of [H] is based on the following 
novel construction of asymptotically free noncommutative random variables which 
emerges as a byproduct of our general theory. (See Corollarv l3.4l below.) Let X and 
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Y be bounded (classical) real random variables. Let X^^^^ and y(^) be independent 
N-hy-N diagonal matrices with diagonal entries that are i.i.d. copies of X and Y, 
respectively. Let H^^^ be an N-hy-N deterministic (!) complex Hadamard matrix. 
Then X^^^ and are asymptotically free. It is especially striking 

here that no "extra" randomness has been added to achieve asymptotic freeness. 

The notion of asymptotic freeness of Haar-distributed unitaries and other types of 
random or deterministic matrices has been extensively developed by many authors 
and in many directions. We just mention the papers [6], [7] and |15| as particularly 
important influences on our work. The reader may consult, say, [TJ Chap. 5], |14j 
or |21j for background and further references. 

To the extent we make progress in this paper we do so by side-stepping issues of 
asymptotic freeness almost entirely. Instead we focus on the notion of asymptotic 
liberation (see ^2.2]i which is far easier to define and manipulate than asymptotic 
freeness. We mention in passing that the operator-theoretic paper [5] helped to 
push us toward a point of view emphasizing operators conjugation by which create 
freeness, and in particular we learned the term "liberation" from that source. 

Another influence on the paper comes from applied mathematics, specifically 
the analysis of high-dimensional data. See for example the paper [19' , which in the 
applied setting makes use of Hadamard matrices randomized both through ran- 
dom choice of block and randomization of signs through multiplication by diagonal 
matrices with i.i.d. diagonal entries of ±1. We use a similar randomization to cre- 
ate arbitrarily large asymptotically liberating families from a single deterministic 
Hadamard matrix. (See Corollarv 13.21 below.) 

The paper |13) considers problems at least superficially similar to those consid- 
ered here, the simplest of which also have the form yl(^) + [/(^)*s(^)C/(^), but 
the focus is not on finding special distributions for J7^^^ to create freeness; rather, 
with t/*^^-' a uniformly distributed permutation matrix, the goal is to calculate the 
limit under various hypotheses on A*^^^ and B^-^\ using the theoretical framework 
of traffics it is ultimately the purpose of [13j to introduce. The problems considered 
typically fall outside the domain of free probability theory proper. A point of view 
encompassing both that of |13| and of this paper would be very interesting to have. 

Here then is a brief outline of the paper. In fj2] we introduce the notion of as- 
ymptotic liberation, we state a technical result (see Proposition 12 . 61 below) spelling 
out the relationship between asymptotic liberation and asymptotic freeness, and 
finally we state our main result (see Theorem 12 . 81 below) . In f|3]we work out a num- 
ber of corollaries to the main result. In fjJl'^G reduce the proof of Theorem 12.81 to 
a combinatorial statement (see Theorem 14.21 below). In fj5]we complete the proof 
of Theorem 12.81 by proving Theorem 14.21 In the proof of the latter the references 
[23( Lemma 3.4] and J24I have an important role. In fj6] we supply a brief review 
of the relevant notions of free probability and a proof of Proposition 12.61 In f}7] we 
prove a concentration result along the lines of that obtained in i3, based not on 
concentration for Haar measure but rather on concentration for the group of signed 
permutation matrices. Finally, in ^we make concluding remarks and mention sev- 
eral more important references the relevance of which is easier to explain after the 
main work of the paper has been carried out. 
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2. Formulation of the main result 

2.1. Matrix notation. Let MatAr (resp., Matkxt) denote the space of iV-by-iV 
(rasp., k-hy-i) matrices with complex entries. We invariably denote the entry of a 
matrix A in the i"* row and j*^ column by A{i, j). For A S Mat^xf, let A* E Mat^xfe 
denote the transpose conjugate of A and let |A] denote the largest singular value 
of A. Let Iat G Mat AT denote the A^-by-A^ identity matrix. 

2.2. Asymptotically liberating sequences. Let / be a finite index set. For each 
positive integer N and index i g / suppose one is given a random unitary matrix 
f/^^^ e Mat A, defined on a probability space depending only on N. We say that 
the sequence of families 

f -\ CO 

is asymptotically liberating if for ii, . . . ,i£ £ I satisfying 

(1) e>2, i\^ii, ig-x^ig, iiT^ii, 
there exists a constant c(zi, . . . , i^) such that 

(2) |Etr (f/f ^Aif/f > . . . C/f I < c(H, ...,^g) {A,} ■ ■ ■ {A,} 

for all positive integers N and constant matrices Ai, . . . ,A£ G Mat a? each of trace 
zero. We emphasize that c(ii, . . . ,ii) is independent of N and Ai, . . . , Ag, but may 
depend on ii, . . . , (and thus, since / is finite, can be chosen to depend on £ alone). 

The interest of the notion of asymptotic liberation stems from its close rela- 
tionship with the concept of asymptotic freeness. We will describe the relationship 
precisely in Proposition 12.61 below. 

Since readers familiar with free probability might find the condition ii ^ ii 
appearing on line ([!]) jarring, we hasten to make the following simple observation. 

f ') °° 

Proposition 2.3. As above, let < | U^^^ I > be asymptotically liberating. 
Then for zi, . . . , G / satisfying 

(3) £ > 2, ?i ^ 12, ■ • ■ , ie-i 7^ (but equality ii — ii not a priori ruled out), 

there exists a constant c{ii, . . . satisfying ([2]) for all positive integers N and 
matrices Ai, . . . , Ag € MatAr each of trace zero. 

Proof. For any ii, . . . , G / and Ai, . . . ^ Ag G MatAr we temporarily write 

^UM, . . . , AAr) = Etr (C^-)A,C^-)* . . . t/f' A,C/f . 

Now suppose that ii, . . . ,ii G / satisfy ^ and Ai, . . . ,Ai are each of trace zero. 
To prove existence of c(ii, . . . ,ig) we may assume that £ > 3 and ig = ii because 
otherwise there is nothing to prove. We then have a relation 

4^.l^{Au ...,Ag) = cl^l^_^ (^AeA, ltr(AMi)I^, A^, . . . , 
+ ltr(AMi)C^..,_,(A2,...,A,_i). 

Thus we can take, say, 

c(«i, . . . ^ 2c(«i, . . . , ie-i) + c(i2, . . . , ii-i), 
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by induction on I. □ 

2.4. Remark on constants. In this paper we always use c, C ot: K (perhaps 
adorned with subscripts or function arguments) to denote constants in estimates. 
The numerical values of these constants vary from context to context. When we 
need to recall some particular constant we will take care to reference the line on 
which it was introduced. 

2.5. The relationship between asymptotic freeness and asymptotically 
liberating sequences. (See §6. II below for a brief review of relevant notions from 
noncommutative and free probability theory. Terms not otherwise defined here are 
recalled there. We work in a setup fairly close to that of [U Chap. 5].) Let / be a 
finite index set. For each i e /, let Ji be a further finite index set. Suppose that 

for each positive integer N and index i £ I one is given a family of 
deterministic matrices in Mat^v such that 



max 



y(JV) 



< oo 



and on a probability space depending only on A'^ a random unitary matrix 
C/f^'' e Matjv- For a C-valued random variable Z, let denote the essen- 

tial supremum of \Z\. Let A^^'' denote the algebra of N-hy-N random matrices A 
which are defined on the probability space where the family is defined 

and satisfy 

IIMIIoo<0«- 

We equip A^'^'' with the state 

0(^)(A) = ^EtrA 
Then we may speak for each i G I oi the joint law 

rf):C({X,,},e,7.)^C 

of 



viewed as a family of noncommutative random variables and of the joint law 

M^^) : C({{X,,},ejJ,e,) ^ C 

of 



{{ 



..(A')^(A')..(A')* 
^ i ij ^ i 



viewed as a triangular array of noncommutative random variables. We make the 
following assumptions. 

(4) For each i G /, the limit — lim t[^^ in moments exists. 



(5) sup max max 



< OO. 



f ^ oo 

(6) The sequence of families < | U^^^ > > is asymptotically liberating. 
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The following technical result explains the relationship between asymptotic libera- 
tion and asymptotic freeness, thus setting the stage for our main result. 

Proposition 2.6. Notation and assumptions are as described above. Then the limit 
fi = limjv->.oo A'*-^-' in moments exists and is tracial. Furthermore, with respect to 
II, the rows of the triangular array {{X^ are free of each other, and for 

each i I the joint law of the row {Xijjjgj. is Ti. 

We give the straightforward proof of the proposition in Sj6] below. 

2.7. Further notation and terminology. A matrix W G Matjv with entries 
= ^i^i,a(j) for some signs ei,...,eN G {±1} and a permutation a G Sn 
will be called a signed permutation matrix. For a C-valued random variable Z let 
= (E|Z|^')^/'' for exponents p G [1, oo). For random vectors X and Y we write 

X to denote equality in distribution. 
Here is the main result of the paper. 

Theorem 2.8. Let I be a finite index set. For each positive integer N and index 
i £ I let there be given a random unitary matrix u\^^ G Mat^y defined on a 
probability space depending only on N . For distinct i,i' £ I put U^^^ = U^^^^*u'^J^^ 
to abbreviate notation. Make the following assumptions. 

(7) For each positive integer N and deterministic signed permutation 
matrix W G Matjv, one has iw*ui^^w\ . = iu}p\ . 

s.t. i^i' s.t. i^i' 

(8) For each positive integer I, one has 
sup max rnax "v/TV ( C/.'.^M (a, /3) 

M=-\ i.i'Gl a. 0=1 \ I 



N=l a, 0=1 

s.t. i^i' 



< OO. 



( ^ C30 

Then the sequence of families \ \ u['^^\ \ is asymptotically liberating. 



We will derive Theorem 12.81 from a combinatorial estimate having nothing a priori 
to do with free probability, namely Theorem 14.21 below. 



2.9. Remarks and notes on the theorem. 



2.9.1. No claim is made that Theorem 12.81 exhausts the class of asymptotically 
liberating sequences. It remains to find other useful general classes of examples. 
In particular it should be possible to use many other groups besides the group of 
signed permutation matrices in formulating variants of hypothesis ([7]). The group 
of signed permutation matrices was emphasized here because of its convenience for 
applications. 

2.9.2. In hypothesis ([7]) we merely require the law of the random vector 
l^^ii^'l ii'ei to exhibit invariance under a certain action of the group of A^-by-iV 

s.t. i^i' 

signed permutation matrices. Dependence of the family {J7j^^^}ig/ is not ruled out. 
This is in contrast to the usual run of free independence results. 
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2.9.3. Nothing is said in Theorem l2.8l about the asymptotic freeness of the famihes 
{U^^^}i^i themselves. The theorem only concerns the possibility of making other 
families of noncommutative random variables asymptotically free by conjugation. 

2.9.4. If each random matrix U^^^ is Haar-distributed in the iV-by-A^ unitary 
(resp., orthogonal) group and each family 1 1//^'' [ is independent, then 

\ oo 

J/-'^-' I > satisfies the hypotheses of Theorem 12.81 as is easy to check. 
J iei 1 



{{ 



Af=l 



Indeed, hypothesis ^ holds by definition of Haar measure and hypothesis (|8]) 
holds because each entry of a Haar-distributed A^-by-A^ unitary (resp., orthogonal) 
random matrix is approximately a centered complex (resp., real) normal random 
variable of variance 1/iV. Via Proposition 12.61 one then can recover relatively easily 
results of the type pioneered by Voiculescu [H] of which the simplest and most 
familiar is that under suitable hypotheses, say, A'^^^ and B^-'^'> are deterministic 
N-hy-N hermitian matrices which are |-]-bounded in N and have empirical distri- 
butions of eigenvalues tending in moments to limits fiA and /j,_b, respectively, the 
empirical distribution of eigenvalues of yl^^^ -I- tends in moments 

to the free additive convolution jiA^ Hb- 

3. Application to free convolution of classical random variables 
We now specialize Theorem 12.81 so as to recover and generalize results of [12]. 



3.1. Hadamard matrices. We pause to review a special class of matrices. One 

says that H e Mat^r is a Hadamard matrix (resp., complex Hadamard matrix) if 
is orthogonal (resp., unitary) and \H{i,j)\ = 1 ioi i, j = 1, . . . , N . The fc-fold 
Kronecker product 



(9) 



1 1 

1 -1 



e Matjv (A^ = 2*^) 



is the simplest example of a Hadamard matrix when A^ is a power of 2. More 
generally, the discrete Fourier transform matrix DFT^^'' € Matjv with entries 



(10) DFT^" > {i, j) = —= exp -27ri 



N \ N 



is an N-hy-N complex Hadamard matrix divided by v A^. 

We now have the following specialization of Theorem 12.81 the proof of which 
relies on the Khinchin inequality. 

Corollary 3.2. Let I be a finite index set. For each positive integer N, let 
Iji^) g Mat AT be a deterministic Hadamard matrix, let W^-^"^ € Mat at be a uni- 
formly distributed random signed permutation matrix and let {F)l^^}i^i be an i.i.d. 
family of uniformly distributed random N-by-N diagonal signed permutation ma- 
trices independent of W^-^"^ . Then 

(11) (Iw^wj u |:^Tywl u 

U ^ \VN ) \ ' VN J^eIJN=l 

is an asymptotically liberating sequence of families of random unitary matrices. 
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Proof. Let / D / be a set such that / \ / consists of two elements, say a and b. Put 
jjiN) ^ ^(jv) andi^f ) = I^,. For* e /U{&},put L/f ^ = i^f )ii;^T4^(A'). Our goal 

r ^ CO 

is to show that < | U^'^'' > > is asymptotically liberating. In any case, it clearly 

satisfies hypothesis ([7]) of Theorem l2.8l Thus we have only to check hypothesis 
To that end, arbitrarily fix positive integers £ and N, indices a,/3 € {1,...,A^} 
and distinct indices € /. As in Theorem 12.81 put uj^^^ — ul^^*ul^\ Our task 

independently of A^, a, /3, i and i'. 



is to bound the quantity VN U-^\a, (3) 

(However, the bound may depend on i.) Now if a € {i, i'}, then there is nothing to 
prove, because the bound is obviously 1. Thus we may assume that € I Li {b}. 
Then the quantity to be bounded can be written 

(12) II (a, /3) 



. N 

Note that the diagonal entries of the diagonal matrix d\^^ d['^'^ are i.i.d. uniformly 
distributed in {±1}. Now the Khinchin inequality (see, e.g., \^\ for a textbook 
treatment) says that for constants ai, . . . , ajsi S C, i.i.d. random variables ei, . . . , eat 
uniformly distributed in {±1} and p G [1, oo), one has 

(13) ||E.=i«»e.|L<^p(Etil«.P)'^' 

where the constant Kp depends only on p, NOT on iV or oi, . . . ,aM- Clearly, the 
Khinchin constant Kg bounds the quantity (fT^ . □ 



3.3. Question. Could the random Hadamard matrices D^^^H^^^ appearing in 

(fTTj) be replaced by deterministic Hadamard matrices ijf^'* in sufficiently general 
position with respect to each other so as still to get asymptotic liberation? More 
generally we wonder to what extent randomness could be reduced in the construc- 
tion described by Corollarv 13.21 

Proposition and Corollarv 13.21 together allow us to perform the additive and 
multiplicative convolution of a pair of classical random variables without introduc- 
ing any extra randomness. While we state the result for two random variables, with 
the evident modifications it holds for any finite number of random variables. 

Corollary 3.4. Let X and Y be bounded real random variables with distributions 
vx and vy , respectively. Let {J^^^-*}^^! be a sequence of N-by-N Hadamard matri- 
ces. Let and {i^^^^j^^i be independent sequences of N-by-N diagonal 
matrices with independent copies of X and Y , respectively, on the diagonal. Let 
A^^^ be the algebra of random N-by-N matrices with essentially bounded complex 
entries defined on the same probability space as A:(^) and and define the 
state 0^^) on A^'^'^ by (l>^-^\A) = 'E-^tvA. Then the noncommutative random vari- 
ables A(^) and are asymptotically free. Consequently the law of 
X(^) + j^H'^N)Y(N)^{N)* converges in moments to 

Vx ffl I'Y 

and the law of -i-iy(^)y(-'^)ij(^)*x(^) converges in moments to 

Vx ^ J^y, 

where ffl and Kl denote, respectively, additive and multiplicative free convolution. 
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Proof. Let VF*-^-* be a uniformly distributed random A^-by-A^ signed permutation 
matrix independent of X^^'' and Y^^\ It is easy to see that the joint law 

,{N) 



(Xi,X2 



of the pair 



equals the joint law of the pair 

To prove the corollary we will focus attention on the latter pair. 

Let {Xi}^^U {Yi}^^ be an independent family of random variables such thatX^ 
(resp., Yi) has the same law as X (resp., Y) for all i. Since we are only proving a 
result about moments, we may assume without loss of generality that X*^^-* (resp., 
y*^^-') is the A^-by-A^ diagonal matrix with diagonal entries Xi, . . . , Xn (resp., 
Yi, . . . , Yn). Then by the strong law of large numbers we have 

(14) lim — tr(xW)''' = EX*^ and lim — tr(r(^))'= = EF'^ 

N^oo N Af-i-oo A^ 

for every positive integer k. We of course also have 



a.s. 



(15) sup 
^=1 



< \\x\ 



< oo and sup 

N=l 



< lirr 



< oo a.s. 



since X and Y are assumed bounded. 

Let C/}^' = ly(^) and U2 = By Corollary 



the following holds. 



(16) 
For ii, 



The sequence {ul^\U2" 'y^^^i is asymptotically liberating. 



€ {1,2} and matrices Ai, . . . ,Ai G Matjv put 



noting that trivially 



(N)* 



AN) 



AN)* 



N 



|V'.„...,,,(Ai,...,A,)|<Pi]...[A,l. 

Let Xi^^ = X(^) and X^^^ = y(^). For ii, . . . , e {1, 2} and positive integers A^ 
we define a random variable 



satisfying the bound 

(17) Id. 
uniform in A^. Clearly we have 

(18) m(^)(X,...XJ=E0W 
for positive integers A^ and ii, . . . ,ig £ {I, 2}. 



( X^'^^ ■ ■ ■ X^" '\ 

ii,--- ,iit y it ' ' ' ' ' it ' 

-\\Y\U' a.s. 



< 



1^1, 



Now by HH), Uni) and ^ the hypotheses (g]), (O and 
are satisfied almost surely. Thus almost surely the rule 



/i(Xij • • • XiJ = lim 0, 



(N) 



for finite sequences zi. 



of Proposition 



e{i,2} 



specifies the joint moments of a state with respect to which X and Y are free 
and which has the laws of the given random variables X and Y as its marginals. 
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Furthermore and crucially, since fi is uniquely determined by its marginals, in fact fi 
is deterministic. Finally, by ([T7| and (fTSl) we can integrate to the limit, concluding 



3.5. Compressions and Discrete Fourier Transform Matrices. Asymptotic 
freeness can be used to address the product of random projections. Consider inde- 
pendent random matrices Pa,Pp € Mat at that project onto uniformly distributed 
subspaces of dimensions [aN] and [/3N], respectively, where [•] denotes the inte- 
ger part. Recall that the empirical distribution function of a Hermitian matrix 
H e Mat AT, which we denote by F^, is 



where Xi{H), . . . , Xn{H) denote the eigenvalues of H. Then the expected empirical 
distribution function of PaPpPa converges in moments to the distribution with 
density given in pO| below (see Section 3.6 of 21 ). Corollarv 13 .61 below shows that 
the same behavior occurs when the subspaces are the spans of appropriate random 
sets of standard basis vectors and columns of the discrete Fourier transform matrix 
DFT^^^. When the first result of this form was presented in [12], the observed 
behavior could not be immediately explained: the random matrix model behaved 
as if the component matrices were asymptotically free, yet they were not thought 
to be so. Theorem 12.81 explains this phenomenon and allows the general statement 
below. 

We next explain in more detail the setup involving the discrete Fourier transform 
matrix DFT*-^^ described in (fTO|) above. Let Ri,R2 G Matjv be two coordinate 
projection matrices, i.e. diagonal matrices with all diagonal entries equal to or 1. 
Then i?2DFT^^^i?i is the matrix that first projects onto a set of indices, performs 
the Fourier transform, and then projects onto another set of indices in the Fourier 
domain. In turn, DFT*-^'' i?2DFT'-^-'i?i is the composition of the projection onto 
the span of a set of standard basis vectors followed by the projection onto the span 
of a set of columns of the discrete Fourier transform matrix. To make this operator 
self-adjoint, we compose is with its adjoint. Thus, the eigenvalues of 



(equivalently: squares of the singular values of _R2DFT*- 'Ri) describe the rela- 
tionship between a subspace spanned by standard basis vectors and a subspace 
spanned by Fourier vectors. These considerations and the trivial observation that 
-\/]VDFT''^' belongs to the class of complex Hadamard matrices motivate the fol- 
lowing result. 

Corollary 3.6. Let {H^^^}'n=i sequence of N-by-N complex Hadamard ma- 
trices. Define independent random variables X and Y by 

( 1 with probability a a v f ^ with probability /3 



that /z^^-* tends in moments to fi. 



□ 



Fh{x) = : KiH) < x} 



i?2DFT(^)i?iDFT(^)*i?2 





In turn, define X^-^^^ and y(^) 
distribution function of 



as in Corollary \3.2[ Then the expected empirical 



(19) 



iV 
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converges in moments to the distribution with density 

(20) (1 - min(a, /3))5o(a:) + max(a + /?-!, 0)5^(x) + ^^'^27r^(f i^^p ^^[A-,A+] W, 
where /[a_,a+] denotes the indicator Junction on [A_,A+] and 



\±=a + P-2aP± v^4a/3(l-a)(l-/3). 
In particular, we niay take 



in the above statement, thus recovering the main result of jl2| . 

Proof. By applying Corollarv 13.21 we only need to determine /ix Kl jiy- This is the 
standard example of multiplicative free convolution and can be found in Section 
3.6 of [2T]. The resulting distribution function has the density claimed. □ 

3.7. Remark. In all the corollaries mentioned above it is not really necessary to 
have matrices i/*-^' which satisfy the definition of a complex Hadamard matrix 
exactly. We stuck to that case only to make all the statements above simple. The 
reader can easily verify that in all the corollaries it suffices merely, say, to have 
€ Matyv unitary for every N and to have maxf^^;^ \H^-'^\i,j)\ be a bounded 
function of N . 

4. Reduction of Theorem 12.81 to a combinatorial estimate 
4.1. Formulation of the combinatorial estimate. 

4.1.1. Further notation. Let \S\ denote the cardinality of a finite set S. Let N denote 
the set of positive integers. For iV e N, let {N) = {1, . . . , N). For A g Mat^x^, let 
[A] 2 = (trAA*)i/2 denote the Hilbert- Schmidt norm of A. 

4.1.2. Functions of x- class. Fix l,N efi. Let F : (N)^ C be a function. We say 
that F belongs to the x- class (the letter x is an allusion to Khinchin) if for all 

ii,...,ii-ie (TV), j£ {N)\{ii,...,ie-i} and Ae {1} 

one has 

F{ii, . . .,ix-i,j,i\, . . .,ie-i) = 0. 

4.1.3. Functions of xx-dass. Fix ^, /V G N. Let F : (iV)^^ C be a function. We 
say that F belongs to the xx-dass if F belongs to the x-class and for all 

ii,---,i2e-2e {N), e {N)\{ii,...,i2e-2} and A e (€) 

one has 

F{ii, . . . ,Z2A-2, J, j,«2A-l, ■ • -,121-2) = F{ii, . . . ,«2A-2,/, j',i2A-l, • • • ,«2£-2)- 

Lemma 15.2.41 below explains the structure of functions of xX"Class in more detail. 

We have the following estimate implying Theorem l2.8l This is the main technical 
result underpinning the whole paper. 



ASYMPTOTICALLY LIBERATING SEQUENCES 



11 



Theorem 4.2. Fix £,N e N. Let F : (N)^^ ^ C be a function of XX- class. Let 
Ai, . . . , Ae G Mat n be matrices each of trace zero and put 

i 

A{i) = Y[ Ax{t2X-i,i2x) for i = (zi, . . . , i^i) £ {Nf. 

A=l 

Then we have 

I 

< a max |i^(i)| TT lAxh 

^ ' A=l 

for a constant Ce > depending only on £, NOT on N , Ai, . . . , Ai or F. 

We emphasize that the Hilbert-Schmidt norm |-|2 intervenes here rather than the 
largest-singular- value norm |-]. We will prove Theorem 14.21 in Sj5] In the course of 
proving the theorem we will highlight close connections with Yin's lemma 
[23l Lemma 3.41 and Whittle's results [24l . 



(21) 



4.3. Derivation of Theorem 12.81 from Theorem 14.21 Fix an integer £ > 2 and 
a sequence of indices ai, . . . , G / such that ai ^ a2, ■ ■ . , ai-i ^ ai and a\. 
For each N and index A = 1 , . . . , ^ define a random unitary matrix 



(TV) 



C/^^]:C/<^^^ ifA>l. 



In turn, define a (deterministic) function i^(^) ; (A^)^^ — > C by the formula 
i^(^)(i) = E n yf\nx^x.i2x) for i = (zi, . . . , *2,) G {Nf. 

A=l 

Arbitrarily fix matrices Ai, . . . ,Ae G Mat at of trace zero. We have to show that 

e 



(22) 



ie(JV)2« 



Etr(t//^Ui • • • V^/^^A,) < c(ai, . . . , a,) n I^^ 



A=l 



where c(ai, . . . , ai) is independent of N and Ai, . . . , A^. By hypothesis ([7]), for each 
i = (ii, . . . , i2e) G (N)"^^, we have 



2£ 



(23) 



(24) 



FW(i) = j ^'"^^(i) for ei, . . . , Civ G {±1} and 



FW(a(*i), . . . ,a(*2,)) = ^'(^ni) for ^ e ^- 



Formula ([23|) implies that F is of x-class and in turn formula ([24|) implies that 
is of xX"Class. By hypothesis (|8]) and the Holder inequality we have 



sup max 

N=i ie(w>2* 



oo. 



The bound 1^ now follows from (gT]) and the evident bound {A}^ < VN {Aj for 
A G Mat AT. The proof of Theorem 12.81 modulo Theorem 14.21 is complete. □ 
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5. Proof of Theorem 14.21 

5.1. Background on set partitions. We are strongly influenced by the paper 
|17| and recommend it as an introduction to the circle of combinatorial ideas being 
exploited here. (But we do not follow the notation of jL7j very closely.) We also 
remark that while set partitions (especially noncrossing ones) have a large role 
to play in the combinatorial side of free probability theory, our use here of set 
partitions seems to be in a different spirit. 

5.1.1. Basic definitions. Let e N be given. A partition of {£) is a family 11 of 
subsets of (£) such that ^ H, Ijn = (£) and for aU Bi, B2 € H, if Bi f] B2 ^ 0, 
then Bi = i?2- Members of a partition are called its blocks. We denote the family 
of partitions of {£) by Part(i?). Given an f-tuple i = (ii, . . . , i^) £ {N)^, let 

n(i) = {{A e{£)\tx^i}\ie {N)} \ e Part(f), 

which we call the set partition generated by the ^-tuple i. For example, we have 

n(l, 2, 3, 3, 2, 2, 2, 4) = {{1}, {2, 5, 6, 7}, {3, 4}, {8}} e Part(8). 

Given Hi , 112 S Part(^) we write Hi < 112 and we say that Hi refines II2 if for every 
Bi G Hi there exists B2 G 112 such that Bi C B2. (In other words, 112 arises from 
Hi by coalescing some blocks of the latter.) In this way Part(^) becomes a partially 
ordered set. Let Og = {{X} \ X G {£)} ^ Part(^), which is the minimal set partition 
(the one with the greatest number of blocks). 

5.1.2. The sum-product formula. Fix £,N e N. Let fi,. ..,fi : {N) C be func- 
tions. Write 

/(i) = /i(u) • • • h(H) for i = (Ji, . . . , I,) e {NY. 

Then one has 

(25) /(i) = n E n n e Part(£). 

ig(Ar)* Beiiie{N)beB 
s.t.n<n(i) 

This is trivial but deserves emphasis. 

5.1.3. Mobius inversion. Consider the function 

C : Part(£) x Part(^) {0, 1}, C(n, 9) ^ { J l^fj^e, 

viewed as a matrix. It is by definition the incidence matrix of the poset Part(£). 
Since C is upper unitriangular, it is invertible and has an upper unitriangular inverse, 
denoted /i : Part(^) x Part(£) — )■ Z, which is by definition the Mobius function of 
the poset Part(^). By definition of /i one has 

(26) EM(0.e) = { J ll;J- forneParK^), 

e<n 

which is a special case of the Mobius inversion formula for the poset Part(£). In 
fact it is the only case needed here. 
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5.1.4. Crude bounds. It is well-known how to evaluate ^ explicitly. One has in 
particular 

^l{Oe,U) = II - 1)! for H G Part(£), 

Ben 

and hence crudely 

(27) max \^{Oi,U)\<f. 
Furthermore, we have 

(28) |Part(^)|</ 

since there can be no more elements of Part(£) than functions {t) {£). 

5.2. Set partitions and the structure of functions of and xx'dasses. 

5.2.1. Special classes of set partitions. For G N put 

Partx(^) = {n e Part(£) | B D {{A} | A e {£)} = 0}, 

Partxx(2^) = {n e Part^(2^) I nn {{2A- 1,2A} I A e (£)}} = 0}. 

In other words, (i) Part^(^) consists of partitions of {£) lacking singleton blocks 
and (ii) Part^^(2£) consists of partitions of {2£) lacking singleton blocks and also 
lacking certain special doubleton blocks, namely {1, 2},. . . ,{2£ — 1, 2£}. 

5.2.2. Structure of functions of x- class. A function / : (N)''^ — > C is of x-class if 
and only if for all i G (N)^ one has /(i) = unless n(i) G Part^(£). A similar if 
more complicated remark holding for functions of the xx-class will be formulated 
below as Lemma [5.2.41 

5.2.3. Clump-equivalence. We introduce some ad hoc and admittedly very ugly ter- 
minology. Fix e,N eN. Given H e Part(2£), let 

Clump(n) = y Be {2£). 

Ben(i)\{{2A-i,2A}|Ae(«)} 

Note that 

n € Partxx(2^) ^ (n G Partx(£) & Clump(n) = {2£)). 

Given i = (ii, . . . ,i2e) G (iV)^^ and j = (ji, . . . , j2«) G {N)'^^ , we say that i and j are 
clump- equivalent and write i ^ j if n(i) — n(j) and i\ — j\ for A G Clump(n(i)). 

Lemma 5.2.4. Fix £,N eN such that N > 6£. Let F : (iV)^^ C be a function 
of the xx-class. Then: (i) For i G {N)"^^ we have F{i) = unless n(i) G Partj^(2£). 
(ii) Furthermore, F is constant on clump-equivalence classes. 

Proof. Statement (i) follows directly from the definition of the x-class. We turn now 
to the proof of statement (ii) which is tedious but only slightly more difficult. Fix 
i — («!,..., 12^) G {N)"^^ arbitrarily and let II — n(i). We have to show that 
F(i) = F{i') for i' G {N)'^^ such that i i'. By statement (i) we may assume 
that n G Part^(2£), and since otherwise there is nothing to prove, we may assume 
n ^ Part;^p^(2£). By symmetry we may assume that 

n = e U {{2s + 1, 2s + 2}, . . . , {2£ - 1, 2£}} 
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for some integer < s < £ and set Q which is empty for s = and otherwise 
belongs to Part^x(2s). Let t — £ ~ s. We have 

{i' e (7V)2^ I i ^ i'} = {{ti,...,i2s,ki,h,...,kt,kt)e {Nf' 

I distinct ki,...,kt E (N) \ {ii, i2s}}- 

Iterated apphcation of the definition of the xx-class shows that F is constant on 
the set above; the only reason for assuming N > 6£ is to provide enough "elbow 
room" for this iteration argument to succeed. We omit further details. □ 

5.3. The Fibonacci- Whittle inequahty. 

5.3.1. Fibonacci random variables. We say that a real random variable (j) has the 
Fibonacci distribution if (/)^ = (/> + 1 almost surely and Ecj) = 0, in which case </> 



necessarily has variance 1. (Equivalently: Pr ^0 = 

V5 



_ l±x/5 
2 



5=F\/5 
10 ■ 



It follows also 



that 
(29) 



1 



and {E^'^} 



T oo 



fc=l 



(0,1,1,2,3,5,8,...). 



The last observation about moments justifies our reference to Fibonacci. We em- 
phasize that while E(f) = 0, we have E^*^ > 1 for A: > 2. 



Lemma 5.3.2. Fix £,N eN. Let 



For i = (ii, 
(30)En( 



■ , *2<?) e (A^)^*^ we have 



,(j)N be i.i.d. Fibonacci random variables. 



A=l 



A=l 



(t'i2X 



if i2\-l ^ «2A, 
i/«2A-l = i2A 



> 



1 */n(i) ePart^x(2£), 
otherwise. 



Proof. The first equality is clear, and it is also clear that the middle quantity is a 
nonnegative integer. Now if n(i) € Vaity^^{2£) and i2A-i = «2A for some AG (£), 
then by definition of Part^^(2i?) there exists A' G {£) distinct from A such that 
«2A G {«2A'-i, *2A'}- The desired inequality now follows from ([29)) . □ 



Lemma 5.3.3. Let A G Mat^r be a matrix. Let p G [2,oo) be an exponent. Let 
4>i, . . . be i.i.d. Fibonacci random variables. Then we have 

N 

^ A{i,j){(t)^(j), - 5ij) 



(31) 



< {Ah 



for a constant Kp depending only on p. 

We refer to this inequality as the Fibonacci- Whittle inequality. 

Proof. Let Xi , . . . , X^r be independent real random variables such that ||2p < oo 
and EXi = for i = 1, . . . ,N. The main result of the short paper [24 , which is 
derived from the Khinchin inequality by a short elegant argument, gives a bound 

1/2 



(32) 



N 

E 

«j'=i 



Aiz, j)iX,X,-EiX,Xj)) 



N 

E 

ij=i 



\Aii,.i)\' 



\^i\\2p II^Jll2p 
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for a constant Cp depending only on p. (We note also that estimates similar in 
form to ((32|) are in common use in RMT. See for example O Lemma 2.7] which is 
more sophisticated but in the same spirit.) Upon specializing this result to the case 
Xi = <f)i for i = 1, . . . , iV, we obtain a bound of the desired form. □ 
We now take a long step toward completing the proof of Theorem 14.21 



Proposition 5.4. Fix £,N £ N along with matrices Ai, . . . ,Ai £ Mat at. (It is not 

necessary at this stage to assume that the matrices are of trace zero.) Put 

e 



A{i) = Y[A^{i2X-i,i2\) /or i= iii,...,i2e) £ (N)^ 



x=i 



We have 



(33) 



5^ \A{i)\<K',l[lA, 



iG(iV)2*s.t. 



A=l 



n(i)GPart^^(2£) 

where Ki is the constant from line (1311) . 

We hasten to point out that this result is not new. It can in principle (albeit with a 
different constant) be deduced as a corollary to Yin's lemma |23i Lemma 3.4]. (And 
for that matter, the Fibonacci- Whittle inequality above could also be deduced as a 
corollary to Yin's lemma.) We omit the derivation of (|33|) by way of Yin's lemma 
because the following derivation seems easier to us (although to be sure, some 
readers may have different preferences). 



Proof. Since Part^p(-(2) — 0, we may assume that £ > 2. Let 0i, , 
Fibonacci random variables. We have 



be 



J2 \m< E i^(i)iEn 

i=(ji,...,i2f)e(JV>2« 



n(i)GPartxx(2£) 



A=l 



< 



n 



A|2 



A=i \ie{N)je{N) 

at the first step by Lemma r5.3.2i at the second step clearly and at the last step by 
Lemma [5.3.31 and the Holder inequality. □ 
One last technical but quite elementary observation is required before the endgame. 

Proposition 5.5. Fix£,N £n, J C (A^) and fi, . . . , fe : (N) ^ C such that 

E /a«-0 forX£{l). 

ie{N) 

Put 

/(i) = /i(*i)---/^(*£) fori^{i,,...,Zi)£{Ny. 

Then we have 



(34) 



ii,....iie{N)\j 

distinct 



< £■ 



•21 



E 



l/(i)|, 
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whe 



{N)Ij = {i = (zi, . . . , i,) €e (N)' I {A} e n(i) ^heJforXe {£)}. 

Proof. Using successively the Mobius inversion formula (j26p , the sum-product for- 
mula ((25)) and the hypothesis /a(*) = 0, we obtain an algebraic identity 



«1, ■•■,«£) 



ii,...,i^G(A^)\J 
distinct 



i6((JV)\J)* 

s.t.n(i)=Of 



^ Mo.,n) = E M(o.,n)n ^ 



nePart(£) ie((Af)\j)* nePart(f) -Ben \ie{N)\j tes 

s.t.n<n(i) 

Eje(Ar)\,/ ribes if ^ is not a singleton, 



E ^(0^' n) n 1 _ , if B = {6} is a singleton. 

nePart(£) sen ^ieJJo\ j i ^ & 

The desired estimate follows easily, via the crude bounds (^7)) and (^5)) . □ 



5.6. Completion of the proof of Theorem 14.21 Let us recall the setup. We fix 
£, TV e N and F : (N)'^'^ ^ C of xx-class. We fix matrices Ai,...,Ai, e Mat at of 
trace zero and put ^(i) = 111=1 ^a(«2A-i, *2a) for i = (ii, . . . , i2i) € (iV)^^. Beyond 
these data we now also fix a set partition 11 € PartT^(2£). By Proposition 15 .41 it will 
be enough to show that 



(35) 



E ^(i)^(i) 



i€{Nf's.t. 

n(i)=n 



< max iFii) 



E 

ie{Nf's.t. 

n(i)ePart^^(2C) 



l^(i) 



Since otherwise there is nothing to prove, we may assume that 
n e Partp(,(2£) \ Part^^(2^). Furthermore, without loss of generality, on account 
of the matching homogeneity in , . . . , of the expressions on either side of (1351) , 
we may assume that N > 6£. Finally, we may assume by symmetry that for some 
integers s > and t > summing to £ and a set Q which is empty for s = and 
otherwise belongs to Part^x(2s) we have 

e U {{2s + 1, 2s -I- 2}, ... , {2£ - 1, 2£}} = H. 



We thus place ourselves once again in the setting of the proof of Lemma 15.2.41 
Then the set of i G (N)^^ such that n(i) = H equals the set of pairs (j, k*^^)) where 
j = (ji,...,j2s) e (A^)^" satisfies n(j) 6, k = {ki,...,kt) S (TV)* has distinct 
entries, k^^' = (fci, /ci, . . . , fct, fct) g (iV)^* and fci, . . . , fc* e (iV) \ {j} where {j} = 
{ji, . . . ,j2s}- Furthermore, by Lemma [5.2.41 there exists for each j € (N)'^^ such 
that n(j) = 6 a unique constant G(j) such that F(j,k) — G'(j) for all i = (j,k(^') 
such that n(i) = U. We now calculate. 
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ie(Ar)2's.t. 
n(i)=n 



E GO) E Aa,k(2): 



s.t. n(j)=e s.t. n(k)=Ot 



< max |F(i)|. E 
s.t.n(i)=n 

s.t.n(j)=e 



(2)N 



ke((w)\{j})' 
s.t. n(k)=04 



< i^* max |F(i)|. E E 

' s.t.n(j)=e 

< max |F(i)|- y \A{i)\. 

^ ' ie(jv)=^*s.t. 

n(i)ePart^^(2£) 

At the third step we apply Proposition 15.51 The other steps are routine. The proof 
of ((35)) is complete, and with it the proof of Theorem 14.21 □ 



6. Proof of Proposition 12.61 

6.1. Brief recollection of notions of free probability. For the reader's conve- 
nience, we review key definitions. For background see, e.g., ^1] Chap. 5], [M] or [21) . 
Let / be a finite index set, and for each i g / let Ji be a further finite index set. 

6.1.1. Noncommutative probability spaces. A noncommutative probability space is a 
pair (A, (j)) consisting of a unital algebra A (invariably having C as scalar field) and 
a linear functional : .4 C called a state satisfying 0(1^) = 1. The state </) is 
called tracial if (j){ab) = 4>{ha) for all a,b & A, in which case by abuse of language 
the space (.4, (f) is also called tracial. 

6.1.2. Unital subalgebras and bracket notation. Given a unital algebra A^ a subspace 
^ -^0 <^ closed under multiplication is called a unital subalgebra. Given a unital 

algebra A and a family of elements {ai}i^j of A, let C{{ai}i^j) denote the unital 
subalgebra of A generated by the family {ai}i^i. 

6.1.3. Noncommutative polynomial algebras. Let C{{'Ki}i^i) denote the noncom- 
mutative polynomial algebra generated by a family {X.i}ig/ of independent non- 
commutative variables X^. By definition the family of all monomials 

(36) X,, •••X„_ (re {0,1,2,...}, zi,...,ve/) 

is a Hamel basis for C({Xi}ig/). Multiplication in C({Xi}i(=/) on the basis acts 
by juxtaposition and the empty monomial is identified with lc({Xi}i(:/) • Given a 
family {a^j.^g/ of elements of a unital algebra A and / = /({X^j.^g/) e C({Xi}ig/), 
the evaluation /({flijig/) G .4 is defined by substituting Xi = m for all i e /. 
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6.1.4. Freeness. Let (A, (p) be a nonconimutativc probability space. Let {Ai}i(=i be 
a family of unital subalgebras of A. One says that the family is free (with 
respect to ^) if the following statement holds for each integer £>2. 

(37) For all zi, . . . , € J such that ii ^ 12, ■ ■ ■ , H-i 7^ H and all ai, . . . , € A 
such that aj € Ai- and 4>{aj) = for j = 1, . . . , one has (j){ai ■ ■ - ai) =0. 

A family {oijig/ of elements of A is called /ree if the family of subalgebras {C{ai)}i^j 
of A is free in the sense just defined. More generally, a triangular array 
{{^ij}jeJi}iei} is said to be free (more precisely, to have rows free of each other) 
if the family of algebras {C({ay is free in the sense just defined. 

6.1.5. Joint laws. Let {A, (p) be a noncommutative probability space and let {ai}i^i 
be a family of elements of A. The joint law ^-l{a,},^, of {ai}i^j is the state on 
C({Xj}jg/) induced by the push-forward rule M{ai}i(=j(/) = 4>{f{{0'i}iei)) for all 
/ e C({Xj}i£j). Recall that in the presence of freeness, the joint law of a family 
{oijig/ of noncommutative random variables is uniquely determined by the laws of 
the individual variables Oi. More generally, if a triangular array {{aij}jeJi}je/ has 
rows freely independent of each other, then the joint law of the array is uniquely 
determined by the joint laws of the rows {aij}jeJi for i G I. 

6.2. Generalities concerning asymptotic freeness. 

6.2.1. Convergence in moments. Given a collection {/i} U of states on 
the unital algebra C({Xi}ig/), we write lim/x^^^ = ii limAr_>oo M^^H/) = m(/) 
for all / e C({Xj}jg/) and we say that {/x^^^} converges in moments to fj,. (Some 
authors on free probability call such convergence weak but we refrain from doing so 
here since the analogue in classical probability really is convergence in moments.) 
Note that if lim/yf_).oo yu^^' = /x and all the /x^^^ are tracial, then so is yu. 

6.2.2. Asymptotic freeness. Let {(^^^\ (/)^^^)}w=i be a sequence of noncommuta- 
tive probability spaces. For each A'' let there be given a family {o-^'}ig/ of ele- 
ments of A^'^^. We say that the sequence {{ai^^}iei}N=i convergent in law if 
the sequence of joint laws l^^^i/^')-^. ^ is convergent in moments. We furthermore say 

{{'^i^^}i6j}w=i i^ asymptotically free if convergent in law and the family {Xi}ig/ 
is free with respect to the limit law on C({Xj}jg/). The preceding definitions can 
be generalized to triangular arrays in evident fashion. 

Lemma 6.2.3. For each positive integer N let 

/.W :C({{Xi,b6jJ.6z)^C 
be a tracial state with marginal states 

rf ) = M^^^lc({x., : C({X,,},e7.) fori€l. 
Assume the following statement holds. 



(38) 



The limit Ti = lim t} in moments exists for each i € I. 
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Assume furthermore that for each integer £ > 2 the following statement holds. 

(39) For all ii^ . . . ,i£ Cz I such that ii ^ 12, 12 7^ *3; ■ ■ ■ , ie-i 7^ ie, 
and for all fi G C({X.,i . . . , G C({Xi,,j}jgj,^), one has 

lim M^^^((/i - M^^H/i)) ■■■{h- m'^H/.))) = 0. 

A'— >oo 

T/ien i/ie ^imif /i — limjv-i-oo Z^*-^-' i?^ moments exists and is tracial. Furthermore, 
with respect to fi, the rows of the triangular array {{Xy j^gj^jig/ are free of each 
other and for each i d I the joint law of each row {X.ij}j^j. is Ti. 

Proof. Consider the following statement for integers £ > 2. 

(40) For all ii, . . . ,ii € I such that ii 7^ 12, • • • , ie^i 7^ ie 

and fi e C({Xji jljgj,^), . . . , ft £ C{{Xi^j}j(zj.^), the sequence 

{t'-^^Kfi ' ' ' fi)} is convergent and moreover converges to 
if lim ^(^) (/j) = for j = 1, . . . , ^. 

In view of hypothesis it is enough to prove that PO]) holds for all ^ > 2. In 

the case ^! = 2, we have 

M^"''(/l/2) = A^(^)((/l - M^"'H/l))(/2 - M^"'^^)) + A^(^n/l)M'"'H/2), 

whence the result by hypotheses (1551) and For £ > 3, statement PO)) holds by 
induction on £ along with hypotheses (|38| and (p9)) . □ 



6.3. Completion of the proof. To prove Proposition 12.61 it will be enough to 
show that the sequence {m'^-'}^^! satisfies the hypotheses of Lemma 16.2.31 At 
any rate, it is clear that each state /i*^^^ is tracial and that hypothesis (l38l) holds 
for the sequence {/i^^^} by hypothesis ([4]) of Proposition 12.61 It remains only to 
check hypothesis (1391) . To that end, fix ii, . . . , and fi, . . . , fg as they appear in 
statement ([Ml) . For X = 1, . . . ,£ put 



4"^^ = .hiMShejJ- ^trA({7;W},eJ,JI^ e An- 

We have 



trA^^'' = and sup max 

N ^ 



< 00, 



the former assertion being clear and the latter following from hypothesis (O of 
Proposition 12.61 We have an algebraic identity 



^/.(^^(/i-M^^H/i)) • • • (/.-M^^H/.))) = Etr(t/f )AWt/f > . . . uipAf^U^P*) 

following immediately from the definitions. Finally, by hypothesis ([6]) of Proposition 
m along with Proposition 12. 3) one has a bound 

Etr(c/f)4"vr*---t/r4"'<h| <c(n,...,^,)ri [[4"^]" . 

Thus ([39| indeed holds. The proof of Proposition 12.61 is complete. □ 



20 



GREG W. ANDERSON* AND BRENDAN FARRELLt 



7. Concentration 

While the asymptotic hberating property of random unitary matrices t/^^^ that 
are invariant under conjugation by a signed permutation matrix ahows us to de- 
termine the hmiting expected eigenvalue distribution functions of matrices of the 
form 

(41) An + U^^^BnU^^> and AnU^^^ BnU'^^^* A*^ 

for suitable sequences of matrices A^, Sjv € Mat^r, the concentration of the empir- 
ical eigenvalue distribution of such matrices is also of interest. Chatterjee obtained 
the first results of this kind when C/'^^^ are Haar distributed [3l |4], and below we 
show that the analogous results hold for the random unitary matrices addressed 
here. In particular, Theorem 17. II is analogous to Theorem 1.1 in [3]. 

Theorem 7.1. Assume that C/(^) e MatAT is a random unitary matrix and that 

WU^^^W* = for any deterministic signed permutation matrix W £ Mat^v- 

Let A, B & MatAT be arbitrary Hermitian matrices and set 

(42) i/+ =A + C/(^'BC/W* 
and 

(43) i/x = AC/(^)B[/W*A*. 

There exists a constant c > 0, independent of A, B and N , such that for all x G 

N 

/ -Nt'^ 
V128(log7V + c) 

for all t > for the same constant c. The analogous two bounds hold as well for 

Hy,. 

Recall that if X and Y are two random variables taking values in the same space 
A", then the total variation distance dTv{^,Y) is defined by 

(45) dTv{X,Y)^ sup |Pr(X e S) -Pr(y e B)|. 

B£X; B is Borcl 

Proof. We construct a random signed transposition matrix in the following way. Let 
J and K be independent and uniformly distributed on {1, . . . , TV}, let {ei, . . . , cat} 
denote the standard basis vectors for and let £1,62 be independent random 
variables taking the values {±1} with equal probability. Define T £ Mat at by 

y ^ f eieje} + Ej5^je»e* J = K 

y eieje^f + e2eKe*j + Y.i^j,i^K ^ie* J ^ K. 

That is, T either does not transpose any entries and multiplies a random entry by 
±1 or does transpose two entries and multiplies them by ±1. 

We let Wr G MatAT denote the random matrix with uniform distribution on the 
set of signed permutation matrices in Mat at. Set 

H'^ = A + WrU^^^W*BWrU^'^^W;, 

so that = H'j^ , which implies that it suffices to prove the theorem for H'j^ rather 
than H+. Now set 



\av{FH+{x)) < 64^28^. Further, 
(44) Pt{\FhA=c) - ^Fh4x)\ >t)< 2exp 



H'l = A + TWrU^^'^W;TBTWrU^'^^W;T, 
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and note that H'l = H'_^. Set D = TWr - W^, so that 

= [D + Wr)U^^\D + W'*)B{D + Wr)U^^^D 

+{D + Wr)u'^^\D + w;)BDu'^^^w; 

+ (£» + Wr)U'^'^^DBWrU'^'^^W; + DU'^'^^DBDU'-^^W;. 



Since D has rank at most 2, each term in the sum above is of rank at most 2. Thus, 
i-ank{H'^ — < 8. Following [3], we apply Lemma 2.2 of [?] to obtain 

(46) ||Fh;-Fh^||oo<^. 

For any fixed a: <E R, we define the function / which maps Wr to Fh'^ {x). By (|46l) . 

\f{W)~FiTW)\ < A 
for any realizations of T, W and f/*^^^; in particular, 

(47) {E{f{W) - f{TW)fY/^ < A 

for all signed permuation matrices W € Mat at. In order to refer to the main theorem 
of we require a bound on drv (Ti ■ ■ ■ Tk, Wr)- Here we apply Theorem 3.1.3 of [T8] 
(which builds on Theorem 1 of [10]) to obtain 

(48) dTv{Ti ■■■Tk, Wr) < ane-^''/^ 

for all A: > 1 for a constant a > 0. (By choosing a large enough we satisfy the 
condition on k stated in [IHl-) Now the claim for follows by using (|46p . (H71) and 
(|48| in Theorem 1.2 of 0. 

The proof for Hx follows in exactly the same way since H'^ — H'^, defined 
analogously to H'^ and -ff", also has rank at most 8. □ 

8. Concluding Remarks 

If An, Bpf £ Matjv are a sequence of self-adjoint matrices whose spectral distri- 
bution converges in any of several ways to a measure, then the empirical spectral 
distribution of the sum and conjugated matrix defined in Theorem 17.11 also con- 
verges in a corresponding way. It remains to provide a bound of some form on 
the distance between the spectral distribution of the resulting random matrix in 
dimension N and the limiting distribution. For other forms of random matrices one 
may use equations for the trace of the resolvent to prove such behavior. This can 
be seen in [11] for Wigner matrices, for example. The authors of [El [20] prove the 
existence of implicit equations as the dimension of the matrices tends to infinity, 
yet for both the Haar case and the case presented here, these equations have, in 
general, not yet been developed for finite dimensions. The exception is [H], where 
the very special case of coordinate projection matrices was addressed. Coordinate 
projections, however, are far simpler than the general case. 

It is possible that using random signed permutation matrices and Hadamard 
matrices offers a new approach to study the convergence of the empirical spectrum 
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to its limiting behavior. In particular, one can embed an element of Sn into Sn+i 
by inserting the index + 1 at a random location. If Fn denotes the empirical 
spectral distribution of the random matrix model described above in dimension 
N, then one can look at ||Ei^Ar — EF/v+i||oo- If this quantity is sufficiently small, 
then a bound on the distance between the empirical distribution and the limiting 
distribution will follow from Theorem 17. II 

Further questions include the following. What other families of random unitary 
matrices are asymptotically liberating? Do there exist families of unitary matrices 
that are not asymptotically liberating, yet yield sums and products that behave as 
if they are free? 

Finally, Corollaries 13.21 and 13.41 suggest the interesting problem of proving a 
"Hadamard analogue" of the strong convergence result of [9l. New tools may have 
to be developed to solve this problem since there is no obvious analogue of the slick 
method used in [5]. 
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